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DEPARTMENT  OF  THE  ARMY  ETL  1110-2-239 

Office  of  the  Chief  of  Engineers 
DAEN-CWE-HD  Washington,  DC  20314 

Engineer  Technical 

Letter  No.  1110-2-239  15  September  1973 

Engineering  and  Design 
NITROGEN  SUPERSATURATION 

1.  Purpose .  The  purpose  of  this  letter  is  to  provide  guidance 
for  one  evaluation  and  identification  of  those  projects  with 
hydraulic  structures  having  the  potential  to  produce  nitrogen 
super  saturation . 

2.  Applicability.  This  letter  applies  to  all  field  operating 
agencies  having  responsibilities  for  the  design  of  Civil  Works 
projects . 

3 .  References . 

a.  ER  1130-2-334 

b.  ER  15-2-11 

4.  BibliooraDhv . 

-  -  ----- 

a.  ER  1110-2-1402 

b.  EM  1110-2-1602 

c.  EM  1110-2-1603 

5 .  Discussion . 

a.  Nitrogen  super saturation  and  associated  fish  mortality 
due  to  gas  bubble  disease  has  occurred  at  Corps  of  Engineers 
projects  on  the  Columbia  River  in  the  North  Pacific  Division 
( NPD)  and  more  recently  at  the  Harry  S.  Truman  project  in  the 
Missouri  River  Division.  Nitrogen  supersaturation  can  result 
at  any  hydraulic  structure  from  entrained  air  introduced  by 
the  spillway-stilling  basin  action.  As  the  flow  is  subjected 
to  hydrostatic  pressure  in  the  stilling  basin,  a  portion  of 
the  entrained  air  is  driven  into  solution  before  it  has  the 
opportunity  to  rise  to  the  surface  and  escape  into  tne  atmo¬ 
sphere.  A  potential  problem  situation  will  exist  if  the 
characteristics  of  the  flow  within  or  downstream  of  the 
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stilling  basin  are  such  that  the  flow  does  not  have  the 


necessary  turbulence 

dissolved  nitrogen.  Flow  conditions  b 
ducive  to  rapid  equilibration  with  the 
shallow,  turbulent  streams.  The  reaera 
characteristics  of  deep,  slow  moving  r: 
reservoirs  are  relatively  small.  Genei 
suffer  from  gas  bubble  disease  so  long 
depths  below  15  feet.  At  those  depths 
nal  gas  pressures  on  fish  are  approxim; 
fish  swim  to  the  surface,  however,  the 
exceeds  the  external  gas  pressure  on  the  fish 
gas  embolism  or  gas  bubble  disease.  The  toler 
to  levels  of  nitrogen  supersatur ation  depends 
of  exposure  and  the  age  and  species  of  the  fis 
dissolved  nitroaen  levels  referenced  to  surfsc 


turbulence  to  degas  or  purge  itself  of  the  excess 
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resulting  in 
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b.  The  phenomenon  of  nitrogen  supersatur ation  below 
hydraulic  structures  is  complex  and  depends  upon  a  number  of 
factors.  Normally  the  problem  of  nitrogen  super  saturation  h 
oeen  associated  with  aerated  flows  plunging  into  deep  sriili 
basins  with  slow  moving  downstream  flow  conditions.  If  tne 
hydraulic  jump  in  the  stilling  basin  is  a  free  jump,  suffi¬ 
cient  turoulence  should  be  present  to  degas  the  flow  so  that 
dissolved  nitrogen  levels  referenced  to  surface  pressure  wil 


as 


not  exceed  110  percent. 
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oxygen  content  measured  or  by  using  a  calibrated  gas  chrom 
graph.  Techniques  to  estimate  the  percentage  of  nitrogen 
super  saturation  below  a  hydraulic  structure  have  been 
developed  by  NPD  and  by  the  U.S.  Bureau  of  Reclamation  (US 
Inclosure  1  gives  a  summary  of  the  development  and  evaluati 
procedure  for  the  MPD  method.  Inclosure  2  gives  a  summary 
the  USBR  method.  The  technique  developed  by  NPD  was  based 
projects  in  the  Columbia  River  Basin.  The  spillways  are  al 
gate-controlled  ogee  crests  and  with  the  exception  of  The 
Dalles,  they  have  similar  stilling  basin  characteristics. 
The  NPD  method  should  be  used  to  evaluate  the  effects  of 
structures  similar  to  those  in  the  Columbia  River  Basin.  T 
coefficients  for  this  technique  are  based  on  these  types  of 
structures.  The  technique  developed  by  the  USBR  is  more 
general  than  the  NPD  technique  and  utilised  data  from  a  vid 
variety  of  hydraulic  structures.  The  USBR  technique  should 
oe  used  to  evaluate  the  effects  of  structures  other  th 
type  found  in  NPD.  Both  techniques  compute  downstream 
gen  concentration  values  by  considering  such  variables  as 
upstream  concentration,  headwater  and  tailwater  elevations, 
head  loss,  angle  of. the  jet,  residence  time  of  the  bu 
and  oressure  conditions  in  the  basin. 
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a.  It  measurements  or  estimates  indicate  that  a  pctentia- 
for  nitrogen  super  saturation  problems  exists,  then  detailed 
model  studies  of  the  project  may  be  necessary  to  develop 
alleviation  measures.  Assistance  in  the  studies  can  oe 
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me  waterways  Experiment  Station 


:ai  assistance  can 
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Interagency  Steering  Committee  on  Reaeration  Research  and 
Committee  on  Water  Quality  (reference  3b; .  Requests  for  • 
services  of  either  of  these  committees  should  oe  coordma: 
tnroucn  HQDA  ( DAEN  — CWE— H)  WASH 
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5.  Action  Required.  Review  ail 
the  procedures  outlined  in  enclosures  1  and  2, 
potential  for  nitrogen  super  saturation  problems  under 
operating  conditions  including  interim  conditions  during 
str uction  . 
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appropriate  portions  of  Survey-Feasibility  Reports,  Design 
Memoranda,  Detailed  Project  Reports,  etc. 

FOR  THE  CHIEF  OF  ENGINEERS: 
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/rJOMER  3. 


Vvliu, 


Chief,  Engineering  Division 
Directorate  of  Civil  Works 
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The  water  parcel  indicated  in  cross-section  by  the  snaced  area  moves 
through  the  stilling  basin,  decelerating  and  increasing  in  height.  It 
extends  laterally  the  full  effective  width,  u  of  the  stilling  basin  as 
illustrated  in  Figure  3  of  the  main  report. 

We  now  make  the  following  assumptions  for  the  water  oarcel  and 
stil 1 i ng  basin : 

1.  For  that  length  of  spillway  that  is  in  operation 
at  a  given  time,  the  discnarge  is  uniform  along  the 

*Taken  from:  "A  Nitrogen  Gas  (N2)  Model  for  the  Lower  Columbia  River, 
"Final  Report,  Water  Resources  Engineers,  Inc.,  under  contract  to 
US  Army  Corps  of  Engineers,  North  Pacific  Division,  January  1971 


Inclosure  1 


to  to 


STL-  1110-2-239 
15  Sep  78 


crest  (this  is  equivalent  to  assuming  that  the 
properties  of  the  water  parcel  are  constant  along 
any  line  parallel  to  the  spillway  crest), 

2,  The  value  is  the  initial  depth  of  the  spill  before 
the  jump,  "It  is  computed  as: 


7  SZga 

o  ^ 


where 

a  ~  discharge  per  foot  along  the  crest 

2  3  total  reservoir  head  asove  the  stilling 
basin  floor. 


3c  The  only  effect  of  the  roller  which  overlies  the  main 
flew  is  to  increase  the  static  pressure  within  t.ne 
water  parcel  by  an  amount  a^y, . 

4.  A  given  mass  .of  air  y,  is  entrained  as  discrete  buobles 
into  the  water  parcslrtat  the  point  s  ~  0  and  remains  -• 
uniformly  distributed  within  the  water  parcel  as  it  - 
passes  througn  the  stilling  basin. 


The  distribution  of  the  mass  of  air  among  the  various 
bubble  sices  remains  unchanged  during  the  water  parcel's 
journey  through  the  stilling  basin. 


The  dissolved  nitrogen  within 
uniformly  distributed. 

Rata  of  nitrogen  dissolution 
governed  by  Field  an  diffusion 


the  water  parcel  is 


drr  in  the  water  sarcel 
a"as : 


where 

A 

X- 


3  the  mass  of  dissolved  nitrogen  in  the  water 
parcel 

3  rata  coefficient. 
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A  =  total  surface  area  of  the  air  bubbles 
contained  in  the  water  parcel, 

C_,  =  effective  saturation  concentration  of 
“  dissolved  nitrogen  in  the  water  parcel,  and 

C  =  actual  concentration  of  dissolved  nitrogen 
in  the  water  parcel. 

\ 

With  these  assumDtions,  we  can  now  define  the  parameters  M}  A }  and  y_ 
in  equation  A-2  as  functions  of  the  location  of  the  water  parcel  in  the 
stilling  basin. 

Assumption  6  allows  us  to  write  the  mass  M  as  the  product  of 
the  concentration  C  and  the  volume  of  the  water  parcel, 

M  =  {-jy&z)C  (A- 3) 

where  w  is  the  effective  width  of  the  stilling  basin,  i.e.,  w  =  (number 
of  gates  open)  x  (width  per  gate). 

The  saturation  concentration  of  a  gas  such  as  N2  or  02  that  is 
only  slightly  soluDle  in  water  is  governed  by  Henry's  Law  which  states 
that  the  equilibrium  or  saturation  concentration  of  the  gas  in  solution 
is  directly  proportional  to  the  pressure  existing  at  the  gas-liquid 
interface.  In  the  water  parcel  the  pressure  ?  at  an  elevation  2  above 
the  stilling  basin  floor  is 

?  =  ?  +  a  +  cl(v-z)  (A-4) 

o  cr  ! 

where  ?  is  the  atmospheric  (or  barometric)  pressure,  and  the  a  parameters 
are  the  densities  of  the  roller  and  main  flow  as  shown  in  Figure  A-l . 

Hence,  the  saturation  concentration  at  any  elevation  2  in  the  parcel  is 
given  as: 
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3<Zt 


l?0  +  ct^  +  ct(t/-a)]<7* 


(A-5) 


when*  C*  is  the  saturation  concentration  under  one  atmosphere  of  pressure. 
In  aquation  A-5,  the  pressure  term  has  units  of  atmospheres  of  pressure. 


rrcm  equation  A-5,  it  is  seen  that  C„  varies  linearly  with 


It  follows 


that  the  average  or  s ff active  saturation  concentration,  C.7  in  the  water 

«w 

parcel  is  the  value  of  C „  ^  at  mid-depth,  or  at  z  -  v/2:  Thus, 

<2  > 


-  r  3 


l?a  +  *  «(y/2)]C* 


v  -•  «  / 


hoping  that  ^  3  5-y  gives  the  final  form  of  c7  as 


i  3  +  3  ' 

*  q  q" 


(A- 7) 


The  total  surface  area  4  of  the  air  bubbles  in  the  water  parcel 
epencs  upon  the  total  mas 5  of  air  entrained  and,  upon  the  bubble  sics 
distribution.  It  is  not  unreasonable  to  expect  that  one  entrained  mass 
of  air  will  be  distributed  among  the  various  bubble  si z*s  in  a  manner 
similar  to  that  shown  below. 


3  *  fraction  of  total  air 

mass  in  the  water  o a real 

with  cubbies  having  a 

mass  less  than  or  aqua! 

to  «. 
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The  volume  7%  of  an  air  bubble  with  mast  *u  can  be 
gas  law: 
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rr\RT 
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where 

m  =  number  of  moles  of  air  in  the  bubble, 

?.  =  universal  gas  constant, 

T  =  absolute  temperature,  and 
P  -  the  total  pressure  in  the  bubble. 

In  ecuation  A-8,  m  can  be  replaced  by  n-J 28.9  where  28. S  is  the  molecular 
weignt  of  air.  The  diameter  dh  and  the  area  of  a  sphere  are  given  by: 


Now,  combining  equations  A-8  and  A-9,  the  following  expression  results 
for  the  surface  area  A.  of  an  air  bubble  with  mass  n.  : 

D  0 


■10) 


Thus,  if  the  total  air  mass* entrained  per  unit  volume  of  water  at- 1  is 
M, ,  the  total  air  bubble  surface  areas  A',  per  unit  volume  of  water  is 

rL 

found  from  the  bubble  size  distribution  and  equation  A- 1 0  as 


or 
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Finally,  to  get  the  total  bubble  surface  area  in  the  water  parcel  it  is 
necessary  to  integrate  equation  A-l 2  over  the  volume  of  the  parcel  wu 6s, 
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/  /  /  A '  dz  sv  dz 


i  A=i 3) 


os  u  3 


Applying  assumptions  4  and  5  and  substituting  for  a’  from  equation  A- 
gi  ves 
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Replacing  ?  with  equation  A-4  and  integrating, 
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If  the  expressions  for  ,v,  tt,,  and  .4  from  equation  .4-3,  4-7  and 

tt 

5  ■-espect'i ve  1  ■/  are  substi tutad  into  tne  rate  expressions  given  in 

aquation  4-2,  there  results 
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{yw6x)^f  =  (wos)  {  [?Q  +  a^D 


+  (ct-o^y]3/3  -  [?o  +  c^Z?  -  a 


r?  + 
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£Lu]^» 
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We  can  now  write  rate  expression  ^  in  terms  of  the  location  in  the 
stilling  basin  by  using  the  relationship 
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where  v  is  the  velocity  of  the  parcel  and  a  is  tne  discnarge  per  unit 
width  of  the  stillinq  basin.  In  addition,  we  define  a  system  parameter 
K ,  whicn  we  will  call  the  entrainment  coefficient,  as 


v  - 


"A 


=  1  ( Dr^:.  « ,/n^a] 
a  ^  28.9/  L  L  A.  d 
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Substituting 
the  concentre 


C-4-V^ 


equation  A-l 8  and  A-l 9  into  A-l 7  gives  the  expression  for 
tion  change  in  the  water  parcel  as 
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The  solution  is  obtained  as  follows, 
the  midDoint  of  the  stilling  basin  y 


Evaluate  the  pressure  terms  at 
to  obtain 
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where 


a. 


O  -r 
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rewriting  equation  A- 21 »  with  these  substitutions  gives 
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mien  has  the  solution 
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Evaluating  equation  A-24  at  r  *  0,  where  C  equals  the  foreoay  concentration 
ioH  it  -  *  :  where  7  equals  the  stilling  basin  concentration  c~a  yields 
-.re  sail  Iway-stl  11 '  ng  basin  mode!  as 
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ESSENTIALS  OF  N-  MODEL-' 

LOWER  COLUMBIA  &  LOWER  SNAKE  RIVERS 


DEFINITION  SKETCH 


RESIDENCE  TIME  »  t_  -  WDL/Q  -  DL/q  ~  L/V 

R  s  L 


TOTAL  HEAD  LOSS  -  h.  -  H-D-h  -  H-(D+v"  /2g) 

Ls  V  L 


ENERGY  LOSS  RATE 


A  /CR 


AVE.  PRESSURE  ■  ?  ■  ?  +  o(D-Y  )  +  a  (D+Y  )  .  a=>  0.0295  atm. /ft. 


°  4 


4  a*  C<  C 

0 


N,  CONCENTRATION  AT  END  OF  STILLING  BASIN,  L 


1/. 

c  -  PC*  -  (?C*-C,,)  exp  (-K  .  vr-  J) 

3*  t  -  L 

q 

1/3  1/3 


IP 


?  +  a  (D+Y  )  -  j  ?  -  cr(D+Y  ) 

4  °j  L  4  0 


1/3 


:< 


L  A? 


"Y/J  In  ^  y  K^g  (1.023)  T  3  water  Temperature. 


?C*-C 


3  aE  .  a  4  b  are  empirically  determined  frcm  observed  data.  They 


are  shown  below: 


MODEL  COEFFICIENTS 


PROJECT 

£ 

a 

b 

Little  Goose 

1.00 

0.09 

2.45 

Lower  Monumental 

1.00 

0.09 

2.45 

Ice  Harbor 

1.00 

0.30 

1 .00 

McNary 

1.00 

1.00 

2. CO 

John  Day 

1.00 

0.20 

2.10 

The  Dalles 

0.50 

0.80 

2.50 

3onneville 

1.00 

1.90 

1.C0 

1 J  Developed  by  Water  Resources  Engineers,  Inc.  for  the  Corps  in  1971. 
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PREDICTION  OF  DISSOLVED  GAS 
AT  HYDRAULIC  STRUCTURES-^7 
by  Perry  L.  Johnson^  and 
Danny  L.  King-2/ 


Introduct ion 


With  the  increased  interest  in  the  effects  of  hydraulic  structures  cn  the 
dissolved  cas  concentration  of  the  flow,  it  becomes  desirable  to  be  able  to 
oredict  how  particular  structures  operating  under  specific  conditions  will 
chance  the  dissolved  oas  concentration . 

At  existing  structures  a  predictive  ability  would  enable  the  facility  operator 
o  select  the  method  of  release  that  would  have  the  most  desirable  effect  cn 
he  dissolved  cas  concentration  of  the  flow.  Prototype  cata  indicate  that  the 
chanoe  in  the  dissolved  gas  concentrat ion  is  dependent  on  the  type  of  structure 
throuah  which  the  flow  passes,  the  macr.it tide  of  the  discharce,  tne  barometric 
pressure,  and  the  water  temperature.  To  establish  an  coeratinc  criteria  for 
each  structure  based  cn  actual  measurement  of  resulting  dissolved  gas  concen¬ 
trations  would  be  a  difficult  task.  A  predictive  abi 1 ity. could  yield  an 
understand inc  of  a  structure's  potential  and  allow  preparation  for  tne  possicl. 
consequences,  even  if  the  structure  had  never  operated. 


Also,  with  a  predictive  ability  designers  would  nave  an  additional  facto*- 
which  could  be  considered  in  structure  selection.  Depend inc  cn  the  situation 
it  is  conceivable  that  the  dissolved  gas  octential  might  even  control  the 
des  ion.  Planners  could  also  use  a  predictive  ability  to  evaluate  the  potent  i 
effects  of  a  sinGle  hydraulic  structure,  or  a  series  of  nydrau-Iic  structures, 
on  a  river. 

Initially,  the  dissolved  gas  concentration  above  the  structure  (both  oxygen 
and  nitrogen)  is  eoual  to  the  concentration  established  by  the  inflowinc 
stream.  The  nitroaen,  being  relatively  inert,  will  maintain  this  concentra¬ 


tion  for  quite  some  time. 


ne  oxycen,  nowever,  especially  in  tne 


of  a  resevoir,  may  be  depleted  from  tne  decaying  of  ora an  it  material.  Thus, 
if  water  is  released  it  may  be  lew  in  dissolved  oxygen  and  yet  may  cor.ce i vac 
be  hich  in  dissolved  nitroaen.  Furthermore,  the  water  may  oe  high  in  bi ocher 
ical  cxycen  demand  (BOD)  which  would  reduce  the  dissolved  oxygen  concentrate 
in  the  stream  below  the  dam.  Therefore,  the  analysis  should  be  able  to 
evaluate  how  effectively  structures  increase  depleted  gas  concentrat ions  as 
well  as  evaluate  whether  supersaturated  conditions  might  be  created. 

Such  predictive  methods  have  been  developed  for  the  spillways  of  the  U.S.  Ar- 
Corps  of  Engineers  dams  on  the  Columbia  River  (1).  Most  of  these  structures 
are  Geometrically  simi'ar.  They  are  low  head,  run-of-the-r i ver  structures, 
with  cate-control  led  ocee  spillways.  The  still! na  basins  are  also  of  simile: 
desion.  This  similarity  enabled  the  development  of  a  predictive  ana "vs is  tn. 
is  quite  satisfactory  for  the  structures  considered.  The  Bureau  of  Reclama¬ 
tion  has  few  structures  that  correspond  to  these  Columbia  River  dams.  In 
General,  Bureau  structures  vary  widely  in  type  and  size.  Thus,  a  much  more 
generalized  predictive  analysis  is  required  for  sicnif leant  application. 

1/  Reprinted  with  the  permission  of  the  authors 

2/  Hydraulic  Engineer,  Bureau  of  Reclamation,  Denver,  Colorado 

3 J  Chief,  Hydraulics  Branch,  Bureau  of  Reclamation,  Denver,  Colorado 
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As  a  basis  for  development  of  the  analysis,  the  following  data  were  collected; 

1.  Reservoir  water  temperature,  dissolved  oxygen  concentration,  and 

dissolved  nitrogen  concentration  at  the  elevation  from  which  the 
water  is  withdrawn 

2.  Discharce  and  a  record  of  which  gates  or  valves  are  coerating  if 

releases  are  being  controlled 

3.  Tai 1 water  elevation,  temceraturs,  and  dissolved  oxygen  and  nitrogen 

concentrations  in  the  tailrace 
>1 .  Local  barometric  pressure 

3,  Rhotocraphs  of  the  structure  coerating  and  dimensioned  drawings  of  the 
structure's  conf iaurat ion 

3y  fall  of  1973  the  monitoring  program  of  the  Bureau's  engineering  and  Reseat: 
Center  had  reached  in  sites  and  had  observed  24  structures  in  ocerat icn . 
Forty-nine  different  operating  conditions  had  been  studied.  In  aedition  the 
Pacific  Northwest  Region  of  the  Bureau  of  Reclamation  has  closely  studi 
Grand  Coulee  Dam  and  made  observations  at  35  other  sites.  The  Upper  Mi 
Region  of  the  Bureau  has  oerfomed  monitoring  at  V e  1  lowtai  1  Afterbay  Oam, 
Combined,  these  data  provided  an  adeouate  base  from,  which  the  predictive 
analyses  could  be  developed. 

Anal  vs : s 


‘he  outcast  of  gas  transfer  is  described  by  the  ecuat 

c(t)  =  Cs  -  ( Cc  -  Cm  e  -<t 


where  C(t)  =  final  dissolved  gas  concentr at  ion 
Cs  =  saturation  concentration 
=  initial  concentration 
<  =  a  constant  of  orccort ional i tv 


is*  M> 
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structure  is  located,  with  daily  fluctuations  that  result  from  atmospheric 
conditions.  The  effects  caused  by  daily  fluctuations  in  atmosoheric  pressure 
are  not  large  but  they  may  be  significant  and  should  be  considered  in  the 
evaluation  of  Cs.  In  this  analysis  measured  barometric  pressures  were 
used  when  available.  If  measured  values  were  not  available  a  standard  atmos¬ 
phere  was  assumed  and  barometric  pressures  were  comp 


:G  accor 


WO  .  Md  V  ■; 


i he  deotn  of  water  over  the  flow  in  which  oas  's  oeina  dissolved  is  general  1 v 
dependent  on  the  depth  of  water  in  the  stilling  basin.  Thus,  variations  in 
the  tai 1 water  elevation  will  have  seme  effect.  Throughout  this  analysis  a 
water  death  eoual  to  two-thirds  of  tne  basin  depth  was  used  to  compute  satura¬ 
tion  concentrations.  It  was  thought  that  initially  the  fairly  compact  jeo 
from  a  spillway  or  outlet  would  penetrate  to  the  floor  of  tne  stilling  basin. 
The  flow  would  then  be  deflected  downstream  and  out  of  the  basin.  As  the  flow 
moved  through  the  basin  it  would  be  diffused  and  its  velocity  reduced. 

This  diffusion  would  be  linear  and  result  in  a  triangular  pattern  with  the 
average  depth  through  the  diffusion  being  two- thirds  of  the  total  basin  ceoth , 
Bubbles  rising  from  the  flow  and  incomplete  'low  penetration  would  tens  to 
r-cuce  this  average  depth,  but  the  two-thirds  deptn  was  considered  repress'": i- 
tive  and  therefore  used  in  the  analysis*  A  major  point  of  support  for  the 
two-tn irds  depth  assumption  is  the  fact  that  later  aopl ’cations  proved  tne 
assumption  reasonable.  If  the  flow  being  studied  coes  not  penetrate  to  tne 
bottom  of  the  pool  the  maximum  depth  of  flow  penetration  may  be  used  in  this 
calculation  in  place  of  the  basin  cepth. 

Evaluation  of  Cs  is  achieved  by  summing  the  barometric  pressure  and  twe-thircs 
of  the  basin  depth  (expressed  in  mm  of  rig)  ana  dividing  this  total  pressure  by' 
standard  atmospheric  pressure  (750  mm  of  Hg)  to  obtain  the  average  absolute 
pressure  on  the  dissolving  bubbles  in  terms  of  atmospheres.  This  average 
absolute  pressure  is  then  multiplied  by  the  dissolved  gas  saturation  concent'': 
tion  at  sea  level,  for  the  desired  water  temperature,  to  obtain  Cs. 


ine  next  parameter  rrom  ecuation  1  to  be  considered  is  the  time,  t.  It 
representative  of  the  length  of  time  that  tne  inflowing  jet  with  entrain* 
is  under  pressure  in  the  stilling  basin  and,  thus,  the  length  of  time  tn< 
is  being  dissolved  in  the  flow.  Consideration  of  time  revealed  two  doss 
limitations  that  could  control  its  value.  First,  it  would  seem  that  civ* 
sufficient  time  the  entrained  air  bubbles  would  rise  out  of  the  flew  and 
the  dissolving  of  gas.  In  some  cases  it  would  seem  that  an  evaluation  of 
bubble  rise  time  could  be  used  to  represent  time.  On  the  Gther  hand,  sit 
tions  might  occur  where  the  flow  with  entrained  air  would  pass  through  tr 
basin  and  be  deflected  to  a  shallow  death  in  a  fairlv  short  time.  There* 


:ne  ac 


!  a 1  lenc 
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of  time  required  for  the  flow  to 

During  this  analysis  the  assumes  ion  wa 
ime  periods  might  be  critical  in  specific  s ituat 


v  iiiur  u 

c  c  f 


lac 8  trial  a- 

these  time  periods  might  be  critical  in  specific  situations.  For  eac 
condition  and  structure  studied,  t  was  evaluated  for  both  limitations, 
smaller  of  the  two  computed  values  was  considered  applicable  to  the  part 
situation  and  was  used  in  the  remainder  of  the  analysis. 


Subbig  rise  time.  -  Evaluation  of  t  based  on  the  bubble  rise  time,  t- 


wouio  be,  if  strictly  pursued,  a  very  complex  computation  which  would  probabl' 
produce  questionable  results.  The  vertical  dimension  of  the  jet  (thickness  e 
jet  that  the  bubble  would  rise  through)  is  never  constant.  The  time,  t,  base 
on  bubble  rise  time,  tj_,  was  evaluated  by  dividing  the  calculated  vertical 
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thickness  of  the  jet  at  the  tailwater  surface  by  the  terminal  rise  velocity 
the  bubble.  3y  trial  and  error,  it  was  determined  that  an  assumed  0.02P-i' 
(0.7-mnO  diameter  bubble  with  a  theoretical  terminal  velocity  of  0.636  ft, 
(0.2  m/s)  yielded  the  most  consistent  results  with  resoect  to  observed  ore. 
conditions.  Also,  when  an  analysis  was  developed  that  predicted  X  (equation 
from  two  dimensionless  parameters,  it  was  ~ound  that  the  0.02°-inch-d i aneter 
bubo  1  e  via* ded  predicted  values  of  !<  that  were  consistent  with  the  predicted 
values  of  X  oasec  on  trip  basin  re  tent  ion  time. 

Basin  retention  time.  -  Comcutaticn  of  the  flow  re tent  ion  time,  to.  in  toe 
basin  is  acccmo :  i sned  by  dtvidina  the  oath  length  of  the  flew  py  the  averace 
flow  velocity  alcnq  the  path.  The  path  lencth  is  generally  controlled  Py 
basin  snaoe.  The  oath  length  is  the  distance  from  the  point  at  which  the  js - 
enters  the  tailwater  oool  to  the  po in' 
directed  toward  the  surface  and,  tner 
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o  cemoute  the  averace  flow  velocity  over  the  oath 
to  obtain  the-  jet  velocity  at  the  tailwater  surface  (or  at  one  start  of  One 
“low  path)  from  toe  previous  analysis  o*  bubble  rise  t’ina,  ~o  determine  the 
averace  flow  velocity,  the  velocity  at  the  end  of  the  oath  must  oe  found . 
This  is  cone  throuch  the  use  of  figure  1  wnich  is  a  summary  of  information 
from  studies  of  jet  diffusion  py  Yevdjevich  (2;  and  Henry  13).  Observation 
velocity  distributions  in  jet  diffusions  indicates  that  half  of  the  maximum 
velocity  would  oe  an  approximation  of  the  jet's  averace  veloci ty  at  the  en 
the  flow  path .  This  averace  velocity  mart  also  be  evaluated  py  jivic'nc 
0  i  sch  arcs ,  T,  o  /  the  channel  cross  sect’  c  n  a  1  area,  A ,  ,«n  ’  on  *cu't  assume 
corns  lete  d  t  f  fus  ion  of  the  jet.  ~~e  lamer  p-f  the  comouteo  vei-oc’  t  ies  should 


used,  since  the  averace  "let  /e'oc’tv  at  the  ere  of  the  oat.n  sou  1  d  te 
n ipner ,  but  rot  lower  than  the  averace  velocity  throuch  the  -u 1 1  cross  s 
■na  /a  ioc "  t  i  es  at  the  beginning  arc  end  o~  t"e  flow  oat.n  are  then  ava^ao 
Chen  this  averace  is  civ’  -ted  into  the  *  *cw  oath  <  enc  th  to  cotain  the  bas 
~  ’  :w  ret  action  ti-e  : t  o  • ,  As  srev :ousiy  stated,  t-e  /a' .a  :f  t  to  se  as 
in  ecu  at  i  on  1  i  s  tne  sma"  a-  gf  t.ne  two  com  s  u  t  ed  va 1  ues  t  ■  z  r  to 


~'-a  final  term  ;n  ecu  at  ion  1  to  be  evaluated  is  '< .  <  is  un  ’ ;  <e  t  re  othe'* 

terms  evaluated  in  that  it  is  not  directly  represent  at i ve  cr  any  stec"*’t 
physical  parameter .  <  is  a  measure  of  tne  ao lilt y  of  a  part icu ’ ar  structu 

operating  under  a  oarticular  condition,  to  dissolve  cas .  it  is  -ecresente 
f  h  ^  ■•onwpa  m ~  a  i  p  n  tre  i  nran:  i n c  cn°  r^ r a  ,4 *  w  n  i  c  r  she  w  e  s 3  r  2 ~  t  r, ^ 
qas- 1  i  qu '  c  inter~ice  is  replenished. 
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severe  i  cectrs.  ’These  date  i  ro  i  rated  that  oisso  Ived  cxycen  cancan Cr  it  :o" 
’/ 2 r v  w i  t  ■/  - ■-» * . -•  d  1 0 °  ./,ph,*h  of  j  ^ a y /a  *  e 1  j t  e r  i o  d  ’  s  “ o  •  / 2 e  r  ’  t  ,'■* e c ^ - 
concentrations  2  re  fairly  c  c  n  s  z  a  n  t  As  s  c  rr  e  0 1  r  e  r  reservoirs  j  ’  s  5  0  I  *  e  ■■ 1  c  • 
citt  were  collected  only  rear  the  surf  ace  and  rot  ac  tne  witrarswa:  e  :ev  2' 
■  0 1 ssol v-ed  nitrocen  arc  o.xycan  esneentrat ions  are  measured  a 
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reservoir  surface  and  the  withdrawals  are  made  from  deeo  in  the  reservoir,  the 
measured  values  of  the  initial  dissolved  qas  concentrations,  Cr,  are  probably 
more  accurate  for  nitrogen  than  oxygen.  Even  though  dissolved  nitrogen  data 
were  used  as  a  base  for  the  analysis,  application  of  the  analysis  for  observed 
prototype  conditions  indicates  that  resulting  dissolved  oxygen  levels  may  also 
be  predicted. 


Figure  2  shows  that  the  value  of  K  is  dependent  on  two  parameters .  The  1 
is  Kv/X,  the  velocity  head,  Hv,  at  the  tai ’water  surface  civiced  by  the 
flow  path  length  X.  Hv/X  is  an  energy  gradient  parameter  for  the  flow; 
relates  the  amount  of  energy  in  the  flow  to  the  path  length  in  the  basin 
which  the  energy  is  dissipated.  The  greater  the  value 


s-f  U 


■no  'T'  r>  r*  d 


i  n* 


turbulent  the  basin  flow  and  the  larger  the  resulting  K  value, 
length  used  corresponds  to  the  value  of  t  selected.  If  t?  is  aopl i cable. 


t,hen  the  value  used  for  X  would  be  the  path  length  used  to  evaluate  to.  ; 
if  is  aopl  icafcie,  the  oath  length  is  adjusted  to  determine  the  effect i-. 
path'lenoth  for  the  time  interval,  that  is,  the  length  cf  t;me  the  bubb'e: 
remain  in  the  jet.  Flow  deceleration  is  assumed  linear  and  the  ratio  C' 
ti_/to  is  multiplied  by  the  total  velocity  drop  to  detpr-rmro  the  v@' 


alone  the  adjusted  path  lencth.  The  average  velocity  a ;c 
path  is  tnen  computed  (initial  velocity  minus  one-half  the  vel 
tiolied  by  ti  to  determine  the  adjusted  path  length. 


p  n  t  np  a  d  ' 
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The  other  parameter  on  which  the  value  cf  K  is  based  is  a  ,*acic  cf  the  snear 
perimeter  of  the  jet  to  the  jet's  cross-sect ional  area  at  the  tailwater 
surface.  This  term  is  a  measure  of  the  jet  compactness  and  shape.  The  shear- 
perimeter  for  a  jet  is  defined  as  the  length  of  the  jet's  perimeter  ove*-  whig: 
a  shearing  action  is  occurring  between  the  jet  and  the  water  of  the  stilling 
basin  cool.  For  a  free  jet  plunging  into  a  pool  the  shear  perimeter  would 
ecual  the  total  perimeter  of  the  jet,  while  for  a  flow  passing  down  a  chute 
spillway  and  into  a  basin  the  shear  perimeter  would  be  the  chute  width  at  the 
tailwater  surface.  Situations  exist  ‘where  the  walls  of  the  stilling  bas;n  ar- 
offset  from  the  jet  entering  the  basin.  If  this  offset  is  small.  Questions 
may  arise  as  to  whether  the  sides  of  the  jet  should  be  included  in  the  snear 
perimeter.  This  is  a  judgment  factor  and  is  probably  best  handled  by  individ¬ 
ual  consideration.  Another  common  structure  that  might  raise  a  similar 
question  would  be  a  hollow  jet  valve  discharging  into  a  pool.  Althoucn  the 
flow  would  have  a  ring-shaped  cross-section,  only  the  outside  perimeter  shod: 
be  included  in  the  evaluation.  In  general,  if  it  appears  that  signifi¬ 
cant  shear  will  occur  along  the  section  of  perimeter  in  question  then  these 
lengths  should  be  included  in  the  analysis. 


the  evaluation  of  K  from  f 
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ysis.  This  empirical  approach  is  mandatory  because  of  the  complexity 
•lows  being  considered.  Very  few  of  the  situations  studied  have  blear 
red  flow  conditions  that  are  well  suited  for  direct  analysis.  Net  cnl 
the  jets  that  leave  the  spillway  chutes,  the  valves,  and  the  cates  oft 
e  complex,  but  the  sti  1 1  inq  basin  bools  are  equally  complex.  Any  anal 
nese  flow  conditions  would  be  quite  involved  and  the  accuracy  would  be 
tionable.  However,  the  coefficients  resulting  from  this  analysis  do  h 
t ional  basis  and  are  reoresantati ve  of  the  various  physical  parameters 
coefficients  can  be  interpreted  to  yield  additional  insight  into  the 
ificance  of  the  various  factors. 
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Althouah  some  entrainment  of  air  is  needed  for  the  dissolved  qas  uptak; 
occur,  the  amount  of  entrained  air  reouired  seems  to  be  quite  small,  i 
of  the  prototype  structures  releases  were  exposed  only  briefly  to  the  ; 
some  of  these  cases  the  water  surfaces  of  the  releases  were  also  re! at 
smooth.  Thus,  it  is  assumed  that  little  air  was  entrained.  This  assur 
was  verified  by  the  small  quantities  of  air  that  were  observed  return i < 
the  tai 'water  surface.  However,  in  some  instances,  the  structures  with 
apparent  air  entrainment  were  among  the  worst  in  creating  supersaturate 

^  r  n  p  j  m  1  Q  ^ 


e  to 

1  *■  ,,  r 
rtl.  S' 

air . 

;  y  9 1  _ 

mot  ion 


I  e  a  c  o  i  i 


7  i  r  at  '  ? 


Included  with  the  example  is  a  drawing  of  the  structure  (figure  3)  an 
graphs  (figure  A)  of  operation.  The  computations  are  described  step 
All  critical  points  and  all  judgments  or  approximations  are  discussed  and  tne 
results  of  the  analysis  are  compared  to  actual  field  findings.  Results  are 
also  included  for  examples  for  which  the  calculations  are  not  shewn ,  far' a- 
t’’ors  between  the  observed  and  calculated  dissolved  gas  concentrations  may  os 
attributed  to  several  factors.  First,  anc  probably  one -of  the  most  important 
is  that  tne  entire  analysis  was  based  cn  average  prototype  data.  Therefore , 
t:~e  structures  will  fit  the  analysis  better  tha n  others  and  seme  structu'-ss 
■will  yield  more  accurate  predicted  results.  A  second  significant  source  o- 
/ar’ at  ion  would  be  errors  in  measuring  the  prototype  dissolved  gas  concentra¬ 
tions  .  rne  chemical  analyses  used  are  not  completely  accurate,  out  even  more 
important,  samples  may  be  collected  ■v-cm  raoions  that  are  not  representative 
of  tne  total  f  low.  Extreme  errors  of  th;s  sort  may  or  may  not  be  cov  i ous„ 
several  cases,  two  or  more  readings  were  available  which  gave  some  add4 tic 
assurance.  Variations  due  to  errors  in  data  collect  ion  nay  be  sms1!  or  tha;. 
may  be  suite  large,  -oo  Meat  ion  of  tne  ana  ys  *s  and  use  of  the  tracts  may 
a' so  -esu  1 1  in  seme  error,  out  tin  is  ^r^gr  should  be  small.  All  r  actors 
cons  '  dared ,  tha  results  are  /ery  encouraging. 

Ex  am? 'e.  -  S'uicawav.  -  'he  ~o'  Towing  icrcrmat 4  cn  is  known : 
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the  stilling  basin  pool  is  short  and  unobstructed.  Because  of  the  changing 
slope  of  the  flow  surface  as  it  enters  the  stilling  basin,  the  angle  of 
penetration  was  approximated  to  be  25*  below  horizontal.  The  basin  depth  of 
22  ft  (5.7  m }  was  computed  for  the  deepest  portion  of  the  pool.  Finally,  th 
flow  path  length,  X,  of  95  ft  (29  m)  is  approximately  the  distance  from  the 
ooint  where  the  iet  would  attain  sionificant  penetration  to  the  end  sill  of 


the  basin. 


T  * 


will  oe  ceflectec  upward,  the  flow  w 


easoned  that  at  the  end  sill  a  large  pc? 

no  loncer  be  under  1 


was  reasoned  tnat  at  tne  end  sill  a  iarce  portion  or  tne  f;ow 

higher  prsssu' 
)  8  c  orn p  1 8 1  s .  T'n  0  s  8  dD o  r o  x  i  m  s 1  i  0  n  s 
are  quite  rough ,  but  attempts  to  refine  tne  evaluations  would  yield  only 
si icnt  improvements  and  would  call  for  and  indicate  unwarranted  accuracy. 


ana 


<  1  c  ’ 


!  7 1  n ; 


or  oases  in  tne  oasin  w 


The  absolute  dissolved  nitrogen  concentration  in  the  reserve i **  is  evaluated 
the  first  step  in  the  analysis.  This  is  accomplished  by  ’-eferring  to  appro; 
ate  standard  tables  and  obtaining  the  nitrogen  saturation  concentration  for 
the  specific  water  temperature  (4. a  *C)  and  multiplying  it  by  the  relative 
reservoir  dissolved  nitroaen  concentration  (104  percent;. 


Cl  =  (1.04)  (2( 


21.5  mc/L 


Next  the  potential  absolute  dissolved  nitrogen  concentration  for  the  still i 
basin  is  computed.  As  stated  before,  it  is  dependent  on  tne  barometric 
pressure,  water  temperature,  and  basin  deotn.  Two-thiros  of  tne  basin  dent: 
is  assumed  as  the  average  depth  over  the  'low  while  the  gas  is  being  d  is  sol’ 
Using  this  approximation  an  average  pressure  on  the  flow  (in  atmospheres )  i : 
computed  and  multiplied  by  the  absolute  dissolved  nitrogen  concentration 
obtained  earlier. 


c  ,  677,.  2/3(,2Z!1304.S/13.»)  (20.7)  .  27.4 

^  /ou 


/! 


This  term  has  been  adjustea  to  reflect  the  barometric  pressure  and,  thus, 
structure's  elevation.  If  the  barometric  pressure  is  unknown,  a  standard 
atmesenere  mav  be  used. 


"  h  c 


'wo  of  the  terms  (Cs  and  Cr)  of  equation  1; 

C(t)  =  Cs  -  (Cs  -  Cj)  e  -Kt 


have  now  been  evaluated. 


The  time,  t,  that  gas  is  being  dissolved, 


is  th 


next  term  of  interest.  The  bubble  rise  time,  ti,  is  evaluated  first. 

To  do  this,  the  vertical  dimension  of  the  jet  at  the  tai ’water  surface  is 
found.  The  23-foot  velocity  head  yields  a  velocity  of  42.5  ft/s  (13.0  m/s). 
The  discharge  is  then  divided  by  the  velocity  to  obtain  a  total  flow  cross 
sectional  area  for  three  gates. 


50/42.5 


qi  £ 

W  •  w 


/  w  -n  L. 


Assumina  ecual  flew  through  each  results  in  a  flew  cross  sectional  area  of 
27.3  f  1 2  (2.5  m2)  for  a  s  inale  gate.  When  equal  flow  conditions  are 
assumed  for  the  gates,  the  analysis  of  each  individual  gate  is  identical  and, 
thus,  the  analysis  'of  the  flow  for  only  one  gate  will  predict  the  performance 
of  the  entire  structure.  If  the  flow  cross  sectional  area  is  then  divided 
by  the  gate  width  (3  ft)  the  flow  depth  is  determined. 


27.3/3  =  3.5  ft  (1.1  m) 
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Since  the  flow  is  not  horizontal  the  flow  depth  must  be  divided  by  the  cosine 
of  the  anqle  of  penetration  to  obtain  the  vertical  dimension  of  the  jet. 

3.5/cos  25*  =  3.5/0.9063  =  3.9  ft  (1.2  m) 


If  this  distance  is  then  divided  by  the  terminal  bubble  velocity,  a 
rise  time,  t,  is  obtained. 


9/0.596  =  5.5  seccne 


"The  length  of  time,  t,  is  also  evaluated  by  consider 
that  the  flow  is  at  an  effective  depth  in  the  basin 


To  do  this  the  curves 


f lours  1  are  used.  First,  the  flow  path  lenqtn,  X,  is  divided  by 
deotn,  rs0. 

X /' 3 0  =  95/3,5  -  27.1 

The  flow  width  (L0)  is  then  divided  by  the  flow  death. 


L0/30  »  3/3.5  =  3.3  • 

riour^  ^  Is  then  referred  to  end  the  ratio  of  the  nax  ’ Turn  velocity,  7—-,. 
within  the  velocity  distribution  at  the  end  of  tne  flow  oath  tc  the  ;  n  ’  1 1  a  ’ 
flow  velocity,  70,  is  obtained. 

VV0  =  9.35 
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head,  Hv,  to  the  appropriate  flow  path  lenqth,  X,  is  Hv/X.  If  the  time 
interval  used  is  based  on  basin  retention  time,  the  basin  flew  path  lenGth 
(evaluated  from  the  basin  geometry)  is  used.  If  the  smaller  time  results  from 
the  consideration  of  the  bubble  rise  time  then  the  flow  path  length  to  be  used 
is  less  then  the  basin  flow  path  length.  For  the  sample  problem  the  time 
based  on  the  basin  retention  time  is  the  smaller  so  the  initially  determined 
path  length  of  95  ft  (29  m)  is  used.  Therefore, 

Hv/X  =  ?p/95  =  0.295 

For  application  of  figure  2,  the  second  parameter  that  must  be  evaluated  : s 
the  ratio  of  the  shear  perimeter  length  of  the  jet  to  the  cross  sectional  area 
of  the  net .  For  this  problem  the  shear  peHneter  is  the  jet  width  plus  the 
jet  height  for  each  side  or 

8  +  3.5  +  3.5  =  15.0  ft  (4.6  m} 


The  cross  sectional  area  has  already  been  found 
Thus  the  ratio  is 


is.  c/; 


,8  =  0.54 


The  value  of  K  is  0.1  from  figure  2.  The  user  will  note  the  possibi 
interpolation  error.  All  the  terms  may  now  oe  substituted  into  ecus,  un  i  sh. 
a  dissolved  nitrogen  concentration  that  is  not  corrected  for  barometric 
oressure  is  obtained. 

Cft)  =  27.4  -  (27.4  -  21.5)  =  23.4  mq/'L 

If  this  is  then  divided  by  the  saturation  concentration,  the  percent  nitrogen 
saturation  is  obtained. 


23.4/20.7  x  ICO  =  113  percent 


"he  coserved  value  for  nitroaen,  No  was  also  113  percent.  To  obtain  a 
predicted  absolute  concentration,  multiply  tne  predicted  percentage  by 
absolute  concentration  adjusted  for  barometric  pressure. 

( 1 . 13 ) ( 577/760 ) ( 20 . 7 )  =  20.8  mo/L  of  N2 

Considering  dissolved  oxvaen,  we  compute: 

O  =  (0.85 ) (12.9)  =  11.0  mc/L 

where  12.9  mg/L  is  the  saturation  concentration  of  oxycen  at  4.4  *C. 


577  +  2/3(22) (2Q4. 8/13. 55) 


lL  (lZ.9,  =  17.1  mr 


t  =  3.3  seconds 
K  =0.1 
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all  cf  which  follow  from  the  nitrooen  calculations  above, 
tion  1: 


Aoplyinq  eoua- 


r  ( -  '1 

=  17.1  -  (17 

Th  s  o s r c 9 

n t  oxycen  sat 

12 ,g/12 

T'pa  d  r  a  7  a 

:  cbssrvsd  v a 

cyqen,  u?  .vas  also  percent 
dissolved  a as  would  oei 


joroximation  of  the  percent  data 

(ICO)  (23.4  +  12 .0 )/ f 20 . 7  *  12.9) 
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"his  considers  nitrogen  and  oxycen,  which 
: h e  total  dissolved  cas. 


together  comprise  over  99  o 


ever a  1  other  examples  were  calculated  with  the  following  results: 
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2.  The  basic  eauation  developed  to  predict  the  resulting  dissolved  gas 
concentrations  is: 

C(t)  =  Cs  -  (Cs  -  Cl)  e  -Kt 

where  C(t)  is  the  dissolved  gas  ccncentrat ion  created  by  the  hydraulic  stri 
ture,  Cj  is  the  dissolved  gas  concentration  in  the  reservoir,  Cs  is  t n 
saturated  dissolved  gas  concentration  at  a  depth  which  is  two-thirds  o 
maximum  basin  depth,  t  is  represent  at i ve  of  the  length  of  time  during  whip: 
gas  is  being  dissolved,  and  K  is  a  constant  that  varies  with  structure  arc 
operating  condition.  A  method  is  developed  for  prediction  of  the  K  value. 
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FIGURE  3  -  SLUICEWAY  IN  EXAMPLE  PROBLEM 
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